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ABSTRACT 

The perturbations such as rotation and PC-current have been believed to be 
greatly affecting the pulsar radio emission and polarization. The two effects have 
not been considered simultaneously in the literature, however, each one of these 
has been considered separately and deduced the picture by simply superposing 
them, but such an approach can lead to spurious results. Hence by considering 
pulsar rotation and PC-current perturbations together instead of one at a time 
we have developed a single particle curvature radiation model, which is expected 
to be much more realistic. By simulating a set of typical pulse profiles we have 
made an attempt to explain most of the observational results on pulsar radio 
emission and polarization. The model predicts that due to the perturbations 
leading side component can become either stronger or weaker than the corre- 
sponding trailing one in any given cone depending on the passage of sight line 
and modulation (nonuniform source distribution). Further we find that the phase 
delay of polarization angle inflection point with respect to the core component 
greatly depends upon the viewing geometry. The correlation between the sign re- 
versal of circular polarization and the polarization angle swing in the case of core 
dominated pulsars become obscure once the perturbations and modulation be- 
come significant. However the correlation that the negative circular polarization 
associates with the increasing polarization angle and vice versa shows up very 
clearly in the case of conal— double pulsars. The 'kinky' type distortions in po- 
larization angle swing could be due to the incoherent superposition of modulated 
emissions in the presence of strong perturbations. 

Subject headings: polarization - pulsars: general - radiation mechanisms: 
non-thermal 
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INTRODUCTION 



Pulsars which are famed by their highly periodic signals are now universally accepted as 
fast rotating and highly magnetized (mainly dipolar) neutron stars. The coherent curvature 
radiation due to the ultra-relativistic plasma streaming out along the open d ipolar magnetic 



field lines is believed to be responsible for the pulsar radio emissio n (e.g.. iSturrockl If 971 
Ruderman fc Sutherlandlll975l ; iMelikidze et al.ll2000t iGil et al.ll2004f ) . The individual pulses 
from pulsars in general are highly random in strength as well as their appearance in longitude 
within the pulse window. However the average profiles resulted from the summation of several 
hundreds of individual pulses have well defined shapes, and they are unique in most of the 
cases. Further pulsars in general show a "S" shaped characteristic polarization position angle 
(PPA) swing which is attribute d to the underlying dipole field geometry of emission region 
( IRadhakrishnan fc Cookdll969l ). 



T he average profiles in gener al ar e made up of many com ponents. iRankinl (119831 11990 



19931 ). iMitra fc Deshpandd (119991 ) and lMitra fc Rankin! (J2002|) have recognized that the pul- 
sar emission beams have nested core-cone structure. But the conal components in gen- 
eral show asymmetry in th eir location with respect to the central core component. Hence 
Lyne fc Manchester! (119881 ) have suggested that the emission is 'patchy'. The components 
often show asymmetry in their strengths between the leading and trailing sides of the pro- 
files. Further, some pulsar s show the polarization angle that deviates from the standard 'S' 
curve ( IXilouris et al.lll998l ). 



Among the several relativistic effects that have been proposed to understand pulsar 
emission and polarization, the effects of rotation such as aberration and retardation (A/R) 
and polar-cap current (PC-current) perturbation are found to be quite important. Due to 
pulsar rotation the relativistic plasma gets corotation velocity component which is in addition 
to the intrinsic velocity along the dipole field lines. Hence the net velocity of plasma will be 
aberrated in the direction of pulsar rotation. Therefore, an inertial observer tend to see the 
plasma trajectory which differs significantly from the associated dipole field lines, and hence 
affecting the pulsar emission and polarization (Blaskiewicz, Cordes and Wasserman 1991, 
hereafter BCW 1991; Dyks 2008; Thomas & Gangadhara 2007; Thomas & Gangadhara 2010; 
Dyks et al. 2010; Thomas et al. 2010; Kumar & Gangadhara 2012a, hereafter KG 2012a; 
Wang et al. 2012). On the other hand, in yet another artificial models where the corotation 
of the pulsar magnetosphere is ignored and considered only the effect of PC-current on the 
underlying dipole field. The field lines will get curvature in the azimuthal direction due to the 
PC-current induced toroidal magnetic field which is in addition to their intrinsic curvature 
in the polar direction. Therefore, the trajectory of field line constrained plasma becomes 
significantly different from the unperturbed case and hence affect the pulsar emission and 
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polarization (Hibschman & Arons 2001; Gangadhara 2005; Kumar & Gangadhara 2012b, 
hereafter KG 2012b). 



By taking into account of rotation, IBCWI ( 1l99ll ) have predicted that the PPA inflection 
point lags the midpoint of the intensity profile by ~ Ar/r^c, where r^c = cP/2n is the 
light cylind e r rad ius, and c is the velocity of light and P is the pulsa r rotation period. 



Later iDyksl ( 120081 ) have confirmed this behavior. However, iKGl ( l2012at ) have shown that 
due to the combined effect of rotation, modulation and viewing geometry the phase lag of 
the PPA inflection point with respect to the centr al core will become significantly different 
from ~ At/tlc- On the other hand, IKGl (j2012bl ) have predicted that the PPA inflection 



point can even lead the central core due to the PC-current-induced perturbation. Note that 



and retardation effects ( 


Gangadhara & Guota 


2001; 


GuDta & Gangadhara 


2003; 


Dvks et al. 


2004; 


Gangadhara 


2005 


)• 



Further, by considering the rotation IBCWI (119911 ) have predicted that the leading side 
intensity components dominate over the trailing ones. This is due to the fact that the 
curvature of source trajectory on l eadin g side becomes larger than t hat on the t railin g 
side. Later iThomas fe Gangadharal (|2007|) have confirmed this effect and lDyks et al.l ( 120101 ). 



Thomas et all (120101 ). IKGl (I2012ah and 



Wang et al.l (120121 ) have reconfirmed this behavior. 



Although, statistically the cases with l eading component stronge r are more common, the 
converse cases are also quite significant (ILyne fc Manchester! Il988f ). By co nsidering the PC- 



current-induced perturbation on the underlying dipole field IKGl ( 12012b! ) have shown that 
the leading side components can get either stronger or weaker than the corresponding ones 
on trailing side dependin g upon the viewing geometry and modulation. However the afore- 



3_up' 

mentioned prediction by IKGl (l2012bf ) was in the absence of strong rotation effect. 



In literature two types of circular polarization have been recognized: 'antisymmetric' and 
'symmetric'. If the circular polarization changes its sense near the center of the pulse profile 
then it recognized as the antisymmetric— type circular polarization whereas if the polarity 
of the circular polarizati on does not change through out the pulse profile then it recognized 
as the symmetric— type ( IRadhakrishnan fc Rankinlll990f ). However, either antisy mmetric or 



symmetric circular polarization can be associated with the individual components (IHan et al. 



1998c lYou fc Hanll2006t lKGll2012aU bl). Earlier only the antisymmetric-type c ircular polar 



ization was thought to b e an intrinsic property of curvature radiation (e.g.. iMichel 



Gil fc Snakowski 



1990aU bl; IRadhakrishnan fc Rankinlll990l ; iGil et al.lll993l ; iGangadhara 



1987 



1997 



20101 ) and the origin of 'symm etric ' — type circular polarization was speculated to be through 
propagation effects. Recently IKGl (|2012aJ) have shown however that in addition to the anti- 
symmetric, the symmetric circular polarization can also be produced within the framework 
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of curvature radiation if the effects of pulsar rotat ion, nonunif o rm pl asma distribution and 
viewing geom etry are take n into accoun t. Recently IWang et al.l (120121 ) have also reconfirmed 
the findings of lKGl (j2012af ). IKGl (l2012bl ) deduced the similar result by considering PC-current 
induced perturbation instead of pulsar rotation. 

Pulsars show a diverse behavior in circular polarization among which its association with 
the PPA swing is quite important for understanding the underlying geometry of emission re- 
gion. In the case of pulsars with antisymmetric circular polarization, iRadhakrishnan fc Rankin 
( 119901 ) have found a strong correlation between the sense reversal of circular polarization and 
the PPA swing: the sign reversal of circular polarizatio n from negative to p ositive is asso- 



ciated with the increasing PPA sweep and vice versa. iGangadharal (120101 ) has c onfirmed 
the correlation and proposed it as a geometric phenomenon. Further KG ( 2012al lbl). from 
their simulation of polarization profiles, showed that the cor relation ex is ts when the rota- 
tion a nd PC-current perturbations are less sig nificant. But Irian et al.l Jl998hJYou fc Hani 
( 120061 ) have noticed that the sense reversal of circular polarization near the center of pulse 
profiles is not correlated with the PPA swing. However, they did find a strong correlation 
between th e sense of the circular polarization and the PPA swing in double— conal pulsars. 
KGl ( l2012a[ ) speculated that such a correlation in the case of double— conal pulsars can arise 
if the modulations are asymmetrically located in the conal rings centered on the magnetic 
axis. 

Among the several relativistic models proposed by taki ng into account of pulsar rotation 
and PC-current p erturbations, only the models proposed by KGl (I2012al ). IWang et al.l (120121 ) 
and iKGl (|2012bJ) can explain full polarization state of the radiation field. Further they 
are much more realistic in the sense that in addition to strong perturbations the effects of 
nonuniform source distribution (modulation) and viewing geo metry have b een i ncorporated 



as an essential ingredients. However in the models proposed by IKGl (j2012al ) and IWang et al. 



( 120121 ) as a special case considered the effect of pulsar rotati on a n d igno red PC-current 
perturbation. On the other hand, in yet another artificial case IKGl ( j2012bl ) considered the 
effect of PC-current on the underlying dipole field by ignoring the corotation of the pulsar 
magnetosphere. Since both the effects are found to be quite dominant in affecting pulsar 
radio emission and polarization, they have to be combined together. 

In this paper both the rotation and PC-current perturbations are considered simultane- 
ously, and analyzed their combined effect on the pulsar emission. If we consider separately 
the rotation and PC-current perturbations, and add up the results, the resulting conclu- 
sions can become erroneous as we show in next section and hence this work is an important 
one. Although pulsar radiation is generated via some coherent processes, we perform the 
modeling of pulsar radio emission in terms of single particle curvature radiation. Note that 
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although to the first order the single particle approximation is not a bad assumption, but in 
reality some factors influencing the coherence processes may favor or oppose the effects that 
we consider in this work. We present the theory of single particle curvature radiation in a 
rotating PC-current perturbed magnetic field and analyze the polarization state in section 2. 
In section 3, we present a set of simulated pulse profiles, and speculate on the polarization 
properties by comparing the observed pulses. In section 4, we give the discussions and in 
section 5 make the conclusion. 



2. CURVATURE RADIATION FROM ROTATING PC-CURRENT 

PERTURBED MAGNETOSPHERE 

Let us consider a stationary Cartesian coordinate system-XYZ with the origin O located 
at neutron star center as an inertial observer frame (see Figured]). Consider an inclined and 
rotating PC-current perturbed magnetic dipole with an inclination angle a with respect to 
the rotation axis Cl, which is taken to be parallel to Z-axis. The velocity v of the relativistic 
source S, which is constrained to move along the rotating PC-current perturbed dipole field 
line f, is given by 

v = Kcb + fl x r , (1) 
where b = B/\B\ and B = B + B±. And B is the unperturb ed dipole fie ld, B± is the PC- 



current- induced field, and r is the position vector of the source (lKGll2012bl ). The parameters 
f2 = QQ is the pulsar angular velocity and k specifies the normalized speed of the source 
with respect to the speed of light c along the associated perturbed field line. 

The first term on the r.h.s. of Equation ([T]) is the velocity of source along the perturbed 
field lines. The second term is the induced velocity due to corotation of the pulsar magne- 
tosphere. Note that due to the PC-current-perturbation, the field lines which lie above the 
magnetic axis, tend to azimuthally twist towards the pulsar rotation, whereas those whic h 



lie below the magnetic axis, twist in the opposite directions (see Figure 1 in IKQ l2012bl ). 
Therefore, the contributions to aberration of the source velocity v by the above two terms 
add up for the negative sight line impact angle a, but they try to cancel each other for the 
positive a. However, since the aberration of source velocity due to the effect of rotation is 
larger than that due the PC-current-perturbation for the current which is of order Goldreich- 
Julian current, the net aberration will be always in the direction of pulsar rotation, and it 
is greater for negative a. 



Since the relativistic emissions are beamed in the direction of velocity v with half- 
opening angle I/7, observer can receive the emissions only from a selected emission re- 
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gion whose boundary satisfies n ■ v = cos(l/7), where v = v/\v\ and the sight line h = 
{sin(, 0, cos^} with ( = a + a. But the e xact a nalyti cal solutions for the emission point 
coordinates (9 , O ) and (9 e , e ) (see IKGl (I2012al ) and iKGl (j2012bl ) for their definition) of 
the beaming region are hard to find once the effects of rotation and PC-current perturbation 
are considered, and hence we seek numerical sol utions. Note that in finding exact values for 
coordinates 9o and 0o, the values o btai n ed from IKGl (j2012bl ) have to be used as initial guess 
values for fa s t conv ergence unlike IKGl ( ]2012ai ) wherein they have used those derived from 
Gangadharal told ). 



By using the parameters r n = 0.1, P = 1 s, source Lorentz factor 7 = 400, the current 
scale factor q = 1, a = 10°, and a = ±5°, we computed #0 and 0o, and presented in Figure [2] 
as function of 0'. Due to the PC-current perturbation alone (represented by dotted line 
curves) , the emission points in O shift to later phases for positive a and to the earlier phases 
for negative a, whereas they are mostly unaffected in 9. On the other hand due to the effect of 
rotation (aberration) alone (represented by dashed line curves), the emission shifts to earlier 
phases in both 9 and 0. However the phase shift of emission points in caused due to rotation 
alone is larger than that due to PC-current perturbation alone. This can be clearly seen in 
the phase shifts of the antisymmetric point of 0o (phase at which O is equal to 0° or 180°) 
indicated by arrows (styled same as 0o). As a result, in the more realistic case of rotating 
PC-current perturbed dipole (represented by thick solid line curves), the emission points in 
9 shift to earlier phases in both the cases of a by the same amount as that in the case of 
rotating dipole, whereas in 0, they shift to the earlier phases by a smaller amount in the 
case of positive a and by larger amount in the case of negative a. Note that, if the emission 
region is modulated (nonuniform source distribution) in the azimuthal direction which we 
show later, then the resulted intensity components will also show aforesaid asymmetric phase 
shift between the positive and negative a cases. 

We also computed 9 Q and O f° r the case of rotating perturbed dipole by assuming 
9o = 9' + 59 TOt + 59 pcc , and 0o = O + ^0rot + £0 P cc, where 9' and 0' o are the coordinates 
of emission points in the nonrotating unperturbed dipole. The changes in the colatitude 
<5#rot = #o iro t — ^0 an d £# P cc = #o,pcc _ 9' and those in the azimuth <50 rot = 0o,rot ~ 0o an d 
^0pcc = 0o,pcc — 0o are due to the aberration and PC-current perturbation, respectively. The 
parameters ^rot and ^o pcc are the coordinates after separately considering the perturbations: 
rotation and PC-current, respectively, and similarly 0o, ro t and 0o, pcc - We find that thus 
obtained 9$ and 0o more or less match with the thick solid line curves (not shown in the 
Figure). Therefore, the phase shifts as well as the changes in the magnitude of emission point 
coordinates #0 and 0o due to the two separate perturbations (rotation and PC-current), 
simply add up when the two effects are combined together. 
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The acceleration of source is given by 
|o| utf |o| a# 

where we have used the expression of arc length of the field line ds = \b\d8 = kc dt, and the 
total derivative dF/d9 = dF/d9 + (dF/d(f>)(d(j)/d9), where F stands for b, n, etc. The first 
term on r.h.s. of Equation (J2]) is the acceleration of bunch due to curvature of the PC-current 
perturbed field lines. Note that it includes both the intrinsic curvature due to the dipolar 
field and the induced curvature due the PC-current perturbation on the dipole field. The 
small change in the source speed due to motion along the perturbed field line is represented 
by the second term and it is tiny among other terms. The third term is the rotationally 
induced acceleration due to the Coriolis force and is in the direction of pulsar rotation. The 
last term is the induced acceleration due to the Centrifugal force which is acting away from 
the rotation axis. The PC-current- induced acceleration which is included in the first term 
becomes much important at higher emission altitude, larger pulsar angular velocity, and 
smaller a. On the other hand, the Coriolis and the centrifugal accelerations with the Coriolis 
term being the dominant one become much important at higher emission altitude, larger 
pulsar angular velocity, but at larger a. 

By using the parameters, which are the same as in Figure (H we computed the radius of 
curvature p^f 3 /|t>xa|asa function of 0', and plotted in upper panels of Figure El The 
PC-current perturbation (represented by the dotted line curves) leads to larger curvature 
with respect the unpertu rbed ones ( thin solid line curve) leaving the symmetric point of p 



to be mostly unaffected (jKGll2012bl ). On the other hand, rotation induces an asymmetry 



into the curvature between the leading and trailing sides of <p' = 0°, since the leading sid e 



trajectory becomes more curved and p maximum shifts to the trailing side (e.g.. lKGll2012al ). 
Since the perturbation due to PC-current in the case of positive a opposes that due to the 
corotation, the net p in the rotating PC-current perturbed dipole (thick solid line curve) 
becomes larger than that due to rotation alone, and vice versa in the case of negative a. 

To assess the possibility of deriving p for the rotating PC-current perturbed dipole by 
considering separately the perturbations due to PC-current and rotation, analogues to the 
#o and 0o in Figure [21 we considered the net curvature as k = k' + Sk rot + 6k pcc , where k' is 
the curvature of the unperturbed dipolar field lines, dk rot = k rot — k' and Sk pcc = k pcc — k' . 
The fc rot and fc pcc are the curvature vectors in the cases of rotating dipole and PC-current 
perturbed dipole, respectively. The resultant p can be derived as 



, , * a, a — a \ n -i/2 

1 1 / (X O ro t O ■ Opcc Q'vat ' *^pc 



-71 + — + — ~ 2 — + 

P Pmt Ppcc V P Prot P Ppcc Prot Ppcc 



(3) 
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where p', p rot , and p pcc are the radii of curvature and of, a rot , and a pcc are the unit acceler- 
ation vectors in the cases of non-rotating dipole, rotating dipole and PC-current perturbed 
dipole, respectively. Thus obtained p, by adding the two separate perturbations computed 
independently, is superposed in p panels of Figure [3] (see the dot-dashed line curves). We can 
see that, it is significantly different from the actual p of the rotating PC-current perturbed 
dipole. 

The polarization position angle ip of the electric field of radiation due to relativistic 
sources, defined as the angle between the radiation electric field and the projected spin 
axis on the plane of the sky, can be computed by knowing the acceleration of the radiation 



source: tan-0 = e ■ a/fen ■ a), where en (pro 



ected spin vector) and e± are unit vectors in the 



directions perpendicular to n ( lGangadharall2010[ ). By using the parameters which are same 
as in upper panels of Figure EJ we computed the position angle ip as a function of 0', and 
plotted in lower panels of Figure [3j The perturbation due to the PC-current (represented 
by dotted line curves) causes the ip curve to shift upward with respect to the standard RVM 
curve (nonrotating dipole, thin solid li ne curve) while its inflection point remain unaffected 
( jHibschman fc Aronsll200ll ; lKGll2012bl ). On the other hand the rotation (aberration) causes 
the ip curve to shift upward or downward, depending upon sign of a, and always shifts the 
inflection point to the trailing side. Therefore in a more realistic case of rotating PC-current 
perturbed dipole the net ip will be shifted upward with respect to the one in rotating case, 
while the inflection point lies mostly at the same phase as in the rotating dipole. 

We have also computed ip for the rotating perturbed dipole by assuming ip = ip'+5ip Iot + 
8ippcc, where ip' is due to the nonrotating dipole, 5ip TOt = ip TO t — ip', and 5ip pcc = ip P cc — 4>' with 
^ rot and ipp CC are being the position angles after considering the perturbation separately due 
to rotation and PC-current, respectively. We find that thus obtained ip in both numerical 
(represented by thick dot-dashed line curves) and analytical perturbatio n theory (represented 
by th in dot-dashed line curves, see the Equations (D13) and (Gil) in iHibschman fc Arons 
20011 ) significantly differ from the actual ip of the rotating PC-current- dipole. Therefore we 
believe that the two effects have to be combined together in deducing the polarization state 
of the radiation field. 



POLARIZATION STATE OF THE RADIATION FIELD 



The radiation emitted by the relativistic accelerating sources will have a broad spectrum, 
and it is given by ( )Jacksonlll998l ): 



1 qe 



iuRo I c 



27T Rq C 



|b| nx[(n-(3)x g| Mt _ ft . r/c) 



KC 



d9. 



(4) 
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Note that the t ime t in th e above equation has to be replaced by the expression given in 
Equation (6) of KG (2 012a ). and the parameters ft by Q t and <j) by the expres sion given in 
Equation (11) of iKGl ( I2012bl ). We solve the integral using the method given in iKGl ( I2012af ) 
and find the polarization state of the radiation field in terms of the Stokes parameters /, Q, 
U and V. 



3.1. Emission from Uniform Distribution of Sources 

By assuming an uniform distribution of sources throughout the emission region, we com- 
puted the radiation field from the beaming region and its polarization state, and presented in 
Figure HI Since the magnitude of rotation of contour patterns of the total intensity and the 
linear polarization in (9, ft)— are more or less the same as that of the circular polarization, 
we present only the contour patterns of circular polarization. However, the total intensity 
and linear polarization will be maximum at the beaming region center and fall considerably 
towards the bo undary wit h a slightly larger emission towards the larger curvature region 



(see Figure 3 in lKGll2012aU bh 



Due to pulsar rotation, the contour patterns of the circular polarization (panels b, b') 
gets rotated in (9, 0)-plane with respect to those in the nonrotating dipole (panels a, a'), and 
the rotation is from the trailing side to the leading side for both the signs of o. Further there 
arises an asymmetry in the strength of the positive and negative polarities of the circular 
polarization in such a way that the negative circular gets quite stronger. On the other hand 
due to the PC-current perturbation, the rotation direction of the contour patterns of the 
circular polarization (panels c, d) becomes opposite between the positive and negative a 
cases. That is for the positive a it is opposite to that due to the effect of rotation, whereas it 
is in the same direction for the negative a. However the magnitude of rotation of the contour 
patterns due to the PC-current is smaller than that due to the pulsar rotation. Also, due to 
the PC-current, the positive polarity of the circular polarization becomes stronger than the 
negative polarity for both ±<r. Therefore, in the case of rotating PC-current perturbation 
(panels d, d'), the rotation direction of the contour pattern of the circular polarization will 
be in the same direction as that in the panels (b) and (&') but with the lower magnitude 
of rotation for the positive a than for negative a. Further, due to the opposite selective 
enhancement of either the positive or negative polarity of the circular polarization by the 
two perturbations (rotation and PC-current), the circular polarization becomes more or 
less have same strength between the positive and negative polarities, in similarity with the 
nonrotating case. 



Emissions from the beaming region due to different plasma bunches will be incoherently 



- 10 - 



added at the observation point if they are separated by a space larger than the radiation 
wavelength. Hence, the resultant emission that the observer receives will be th e sum of the 



i ntens ities from the different plasma bunches. Using the expressions given in iGangadhara 



( 120101 . from Equation 33 to Equation 36), we computed the polarization state of the emitted 
radiation due to uniform distribution of sources, and they are presented in Figure [5j With 
the uniform distribution of sources, the PC-current perturbation alone (represented by the 
dotted line curves) does not affect the symmetry of the total intensity / between the leading 
and trailing sides of = 0° which is similar to the nonrotating dipole, whereas the rotation 
introduces an asymmetry (represented by the dashed line curves) with the leading side 
becomes stronger due to an induced larger curvature of source trajectory on the leading side 
(see Figure [3]). In the more realistic rotating perturbed dipole, there remains an asymmetry 
similar to the case of rotating dipole, however with the emission that is significantly differs 
from that in the cases of the rotating dipole and the PC-current perturbed dipole when 
considered separately. 

Due to an incoherent addition of emissions from the different bunches within the beam- 
ing region the magnitude of linear polarization L become bit smaller than that of total 
intensity /. However, the profile of L more or less matches with the corresponding /. The 
net survived circular polarization V in the case of rotating perturbed dipole (represented by 
the thick solid line curves) becomes very tiny because of the two perturbations (rotation and 
PC-current), which selectively enhance the opposite polarities of V as shown in Figure HI 
The position angle is increasing in the case of positive a whereas it is decreasing in the case 
of negative a with the inflection point always shifted to trailing side. 



3.2. Emissions Due to Nonuniform distribution of Sources 



We considered a Gaussian modulation function that is given in IkgI (boi2aU rl) to model 
the nonuniform distribution of sources in the emission region. We find the polarization state 
of the radiation field f rom the nonuniform distribution of sour ces using the expressions given 
in 



Gangadharal (120101 ) (see equations 38 of iGangadharal 120101 ) . 



3.2.1. Emission with Azimuthal Modulation 

By considering a modulation in the magnetic azimuth with the peak at the meridional 
plane, we simulated the polarization profiles affected by the rotation and PC-current per- 
turbations, and plotted in Figure M along with the profiles of unperturbed emissions. We 
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chose r n = 0.05, f$ — 1, <pp — 0°, cr^ = 0.1 and the rest parameters the same as in Figure [5j 
We can see that due to perturbations both the emission and polarization significantly get 
affected in phase and magnitude, however, the maximum of L mostly remain unaffected. 
Pulsar rotation causes the intensity components to get shifted to the earlier phases with 
respect to the fiducial plane, whereas the PPA inflection points to later phases for both the 
signs of a. On the other hand PC-current causes the intensity components shift to later 
phases and the PPA inflection point to earlier phases for the positive a and vise versa for 
the negative a. 

The net phase shift of intensity component and that of the PPA inflection point after 
combining the two perturbations are found to be —2°. 16 and 7°. 95, respectively in the case of 
a = 5°, whereas they are found to be —6°. 69 and 8°. 90, respectively in the case of a = —5°. 
On the other hand, the sum of the phase shifts of the intensity components caused separately 
by the rotation and PC-current are found to be —2°. 10 and —6°. 87 respectively in the cases 
of a = ±5°, whereas that of the PPA inflection point are found to be 6°. 44 and 9°. 59, 
respectively. Therefore, the absolute relative difference between the phase shift of intensity 
component that resulted when the two perturbations taken together and that due to the 
sum of two separate perturbations, is about < 3% in both the cases of a, whereas that for 
the PPA inflection point is found to be about 19% in the case of a = 5°, and 8% in the case 
of a = —5°. Note that the net phase shift of the intensity component becomes slightly larger 
than that of the net modulation / (which is about — 1°.92 and —5°. 71, respectively in the 
cases of a = ±5°) due to an induced asymmetry in the net radius of curvature about the 
peak location of modulation. 

Although circular polarization of opposite polarities from the background unmodulated 
emission roughly cancels out, a net circular however with an asymmetry between the opposite 
polarities survives in the presence of modulation. This is because of an asymmetry in the 
magnitude of rotation of the emission pattern with respect to rotation phase in such a way 
that larger ro tation magn itude towards the inner rotation phases as compared to that on 



outer phases (IKG l2012al ). Hence, it results in the selective enhancement of the leading 
side circular over the trailing side circular. However, the asymmetry between the opposite 
polarities of the circular polarization becomes smaller for the positive a than that for the 
negative a. This is due to an opposite behavior of PC-current in introducing an asymmetry 
between the opposite polarities of V between the ±<r wherein it selectively enhances the 
trailing positive circular for +a whereas the leading positive circular for —a. On the other 
hand pulsar rotation causes the selective enhancement of the leading polarity of V, i.e., 
negative circular for +a and positive circular for —a. 

To assess whether we can estimate the net phase shift of intensity components by adding 
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the phase shift due to each perturbation computed separately, we chose emission altitude 
r n = 0.1, and two cases of modulation: narrower (o^ = 0.1) and broader (cx^ = 0.3). 
The simulated pulses are presented in Figure The net phase shift of the component 
after combining the two perturbations in the cases of cr^ = 0.1 and 0.3 are found to be 
—4°. 40 and —8°. 33, respectively, whereas the PPA inflection point phase shifts are found 
to be 15°. 70 and 16°. 23, respectively. On the other hand the sum of the phase shifts of 
the intensity components caused due to the two perturbations considered separately are 
found to be —4°. 09 and —5°. 77 respectively wherein the relative difference are about 7% and 
30% respectively. For PPA inflection point the corresponding phase shifts are found to be 
14°. 52 and 16°. 52, respectively, wherein the absolute relative differences are about 7.5% and 
2%, respectively. Hence the actual phase shifts of intensity components and PPA inflection 
point will be significantly different from those obtained when the two perturbations dealt 
separately. 

In the case of = 0.1, the asymmetry in the strength of circular polarization between 
the opposite polarities become larger as compared to that in the lower altitude emission (see 
Figure E]). This is due to an increased asymmetry in the magnitude of rotation of beaming 
region emission pattern between the inner and outer phases wherein the small amount of 
rotation occurs at the outer phases. In the case of = 0.3, only the leading negative 
polarity of the circular polarization survives due to aforementioned reasons but this time 
it is enhanced more due to broader modulation and hence broader pulse width. Hence it 
results in the "symmetric" -type circular polarization. 

As a case of modulation which lies symmetrically on either side of the meridional plane, 
we choose the modulation peaks located at <pp = ±30° for a = 5° and <pp = 180° ± 15° for 
a = —5°, and the simulated profiles are given in Figure [HI Note that only the combined 
case of rotation and PC-current perturbations is given in all panels. Although an observer 
encounters mostly the same plasma density (modulation strength) between the leading and 
trailing sides, the net modulated intensity component on the trailing side becomes weaker 
than that on the leading side due to the induced larger curvature of source trajectory on the 
leading side. Further the leading side component becomes broader than the trailing one due 
to an induced asymmetry in the sight line encountered modulation and that in the gradient 
of radius of curvature. The behaviors of linear and circular polarization are same as in the 
Figure [6] except the enhancement of trailing part of the circular polarization in the trailing 
side components. 
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3.2.2. Emission with Polar Modulation 

By considering a modulation with plasma density gradient in the polar direction we 
present the effects of rotation and PC-current perturbation on pulsar emission and polar- 
ization. In Figure |9l we present the simulation of hallow cone emissions surrounding the 
magnetic axis with modulation peak at 9p = 2°. Since the minimum of 9q is about ~ (2/3) a 
which is about 3°. 3 for o = ±5°, observer's sight line cuts the hollow cone emission only 
once in each pulsar rotation, and hence it results in a single component profile. For the sake 
of comparison, the cases of before and after consideration of the perturbations are all shown 
in the Figure M, and the parameters normalization and the line representations are the same 
as in Figure Note that, unlike the case of azimuthally modulated emissions (see Figure 
[6]), the maximum of modulation strength that the observer encounters in the lab frame is 
not same in all the cases as the minimum of 9$ slightly get affected due to the perturbations. 
But, similar to the case of Figure [6], the emissions get affected significantly before and after 
considering the perturbations due to the induced differences in p at the peak locations of 
/. The net phase shifts of a component after combining the two perturbations in the cases 
of o = ±5° are found to be —5°. 79 and —6°. 29, respectively. On the other hand the net 
phase shift obtained by adding the phase shifts due to the rotation and PC-current when 
considered separately in the cases of a = ±5° are found to be —5°. 32 and —6°. 83, respec- 
tively. Hence the relative differences are about 8% and 8.5%, respectively with respect to 
the combined case of the perturbations. Note that this relative differences in the intensity 
phase shifts become larger at higher altitude (~ r n = 0.1). 

The maximum of normalized linear polarization L is more or less the same in all the 
cases, similar to Figure M Although the position angle significantly gets affected after com- 
bining the two perturbations, its inflection point lies at roughly the same phase as that due 
to the rotation alone. The circular polarization V becomes symmetric type with the opposite 
polarities survive due to an opposite direction of rotation of emission pattern in the beam- 
ing region. The net circular polarization after combining the two perturbations will also be 
symmetric type with the survival of polarity as that due to the pulsar rotation alone. 

By considering a hallow cone modulation with peak at 6p = 3°. 6, we have analyzed a 
case where the sight line cuts the hallow cone emission twice (see Figure [TO]) . Similar to 
Figure El the trailing side intensity component become weaker as well as narrower than that 
on the leading side, due to the induced asymmetry in the curvature of source trajectories 
between the two sides. Since the sight line crosses the central maximum of the hallow cone 
on both the leading and trailing sides, there is a change over of selective enhancement of 
emission over the part of the beaming region with smaller values of 9 to that over the larger 
values 9 on the leading side, and vice versa on the trailing side. Hence, there results an 
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antisymmetric circular polarization with the sign reversal from positive to negative over in 
leading side component and vice versa on the trailing side for the case of positive a. On the 
other hand, it is vice versa for the case of negative a due to the opposite direction of rotation 
of the emission pattern of the beaming region. 

3.2.3. Emission with Modulation in both Azimuthal and Polar directions 

The radiation sources may be nonuniformly distributed in both the polar and azimuthal 
directions in the pulsar magnetosphere. The extreme cases, wherein the modulation is pre- 
dominant in azimuthal or polar directions are already discussed in Sections 3.2.1. and 3.2.2., 
respectively. In this section we present a few more cases where the modulation exists in both 
the polar and azimuthal directions. In Figure [TTJ we considered the cases oq = 0.01 and 
o$ = 0.1 wherein the modulation is effectively dominating in the azimuthal direction over 
that in the polar direction, and o~q = 0.001 and = 0.5 wherein the modulation becomes 
effective in 8. Note that in the case of o~g = 0.01 and = 0.1, even though o~q < a^, the 
modulation becomes effective in <p coordinate than in 6 due to the much larger coverage of <fi 
compared to that of 9, see for e.g., Figures [3] and [U In the case of o~q = 0.01 and = 0.1, the 
emission and polarization properties are similar to the case of a = 10° and a = 5° of Figure 
If)] with an antisymmetric circular polarization wherein the sign reversal is from the negative 
polarity to the positive. On the other hand even though the single modulation is considered 
in the case of o~q = 0.001 and = 0.5 but due to the viewing geometry, much elongation of 
the modulation in the azimuthal direction and squeezing into a narrow cone, the modulation 
encountered by the sight line results in a blended two components like structure. It further 
results into a two components like structure in intensity profile, however, with much weaker 
trailing side due to the larger p. Although circular polarization is still antisymmetric type, 
the sign reversal becomes opposite to the case of oq = 0.01 and = 0.1, i.e., from positive 
to negative, which is similar to the case of a = 10° and a = 5° in Figure [9j Hence in 
the presence of perturbations, it is hard to see the correlation between the sign reversal of 
circular polarization and the PPA swing as both the types of circular sign reversal seems to 
be associated with the increasing PPA swing. Note that, the inflection point of PPA swing is 
not derived in the case of o~q = 0.001 and = 0.5 due to the difficulty in finding it because 
of the kinky nature in the PPA swing. 

By using a = —5°, 0p = 180° and the rest parameters the same as in Figure [UJ we 
computed the polarization profiles and plotted in Figure [121 In both cases of modulation, the 
profiles are similar to those in Figure [11] except with the opposite sign reversal of the circular 
polarization in the respective cases and decreasing PPA swing. Also, similar to Figure [11] due 
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to perturbations no correlation is found between the sign reversal of the circular polarization 
and the PPA swing, as both the sign reversal of the circular polarization, i.e., from positive 
to negative or vice versa, are associated with the decreasing PPA. 

In Figure [13j we presented the simulations for different cases of two Gaussian modula- 
tions symmetrically located in a given ring centered on the magnetic axis. Although, the 
asymmetry in p between the leading and trailing sides is same in all the cases with much 
larger curvature on the leading side, there is a diverse asymmetry in the strength of the inten- 
sity component between the leading and trailing sides. This is due to an induced asymmetry 
in the strength of modulation that the inertial observer encounters between the leading and 
trailing sides. In the case of Op = 4° and cf)p = ±15°, observer encounters a much weaker 
modulation on the leading side than that on the trailing side which overcomes the influence 
of asymmetric p. Hence, it results in a stronger trailing side component than that on the 
leading side. Whereas in the case 9p = 4° and <pp = ±25°, observer encountered asymmetry 
in the modulation strength between the leading and trailing sides becomes smaller compared 
to that in the case 9p = 4° and <pp = ±15°. Hence the leading side intensity component 
becomes stronger than the trailing one as the influence of asymmetry in p becomes much im- 
portant than that in /. In case of Op = 4° and <pp = ±35°, the asymmetry in the modulation 
strength between the leading and trailing sides becomes even smaller, and hence resulted in 
mostly in leading side component. In case of Op = 4° and <pp = ±45°, observer encounters 
higher modulation strength on the leading side than that on trailing side, and hence resulted 
in a single leading side component or a partial cone. 

In all the cases, circular polarization is found to be symmetric type over both the leading 
and trailing side components due the selective enhancement of negative polarity caused by 
the rotation and PC-current perturbation. Note that for example in the case of Op = 4° 
and <pp = ±15°, the emissions from the region, which is closer to the magnetic axis, are 
selectively enhanced, and hence a selective enhancement of inner circular occurs (see for e.g., 
the case a = 10° and a = 5° in Figure [TUl) . The kinks are introduced into the PPA swing due 
to the combined effect of modulation and perturbation caused by rotation and PC-current. 

Similar to Figure [T3"| we selected Gaussian modulations which are symmetrically located 
on either sides of the magnetic meridian plane in a given cone, and the simulated profiles are 
shown in Figure [TH for a = —5°. Similar to the cases in Figure [131 the modulation strength 
that the observer encounters can be quite different between the leading and trailing sides 
in addition to an asymmetric p. Hence it leads to a larger asymmetry in the strength of 
components between the two sides. In the case <pp = 180° ± 10° observer misses out the 
trailing side modulation and hence only the leading side component shows up, whereas in 
the case <f)p = 180° ± 40° it becomes vice versa. For third column, we used the combined 
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modulation that are used in the first two columns (four Gaussian modulations). Even though 
observer encounters a weaker modulation on the leading side, due to larger curvature on the 
leading side the component on the leading side becomes stronger than the corresponding 
one on trailing side. But in the case of last column panels, observer encounters even weaker 
modulation on the leading side than on the trailing side as compared to the case of third 
column. Hence the resultant leading side component becomes weaker than the corresponding 
trailing side one. In all the cases, similar to Figure [131 due to the selective enhancement, 
the circular polarization becomes positive symmetric type over both the leading and trailing 
sides. 



4. DISCUSSION 



The poloi dal PC-current perturbs the underly ing; dipole field by inducing a toroidal 



magnetic field ( jHibschman fc Aronsll200ll ; lKGll2012bl ) whereas the pulsar rotation causes the 
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with modulation can shift the intensity to trailing side and the PPA inflec tion point to the 
leading side for the positive a, and vice versa for negative a (lKGll2012bl ). The influences 
of the rotation and PC-current adds up for the negative a but cancel each other when o is 
positive. Note that the opposite behavior of rotation and PC-current in the case of positive 
a is due to the opposite directions of induced curvature caused by the two perturbations. 

Although the effects of rotation and PC-current perturbation can be understand qual- 
itatively when they are combined together, their quantitative estimate becomes impossible 
when the two effects are considered separately as shown in sections 2 and 3. However, since 
the influence of rotation is mostly larger than that of the PC-current, in the combined cases 
mostly the effects of rotation prevail with quite different magnitudes. For example, the phase 
delay between the PPA inflection point and the central component of pulse profile becomes 
much smaller for the positive a as compared to that for the negative a (see Figure [6]). Hence 
the above said phase delays greatly depend on the geometric parameters which has not been 
reported in the literature. 
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Even though the curvature of source trajectory on the leading side becomes larger than 
that on the trailing side due to the perturbation caused by the rotation and PC-current 
considered together, the net modulated intensity can become stronger either on the leading 
side or on the trailing side (see Figures [T3] and [TJJ depending upon the modulation location 
and the viewing geometry. Statistically, although there is an observational support for the 
leading si de component dominating over the trailing one, the other contrary cases are also 
reported (ILyne fc Manchesterlll988l ). Hence our model provides a much plau sible explana tion 
for the usual asymmetry between the leading and trailing side components. iKGl (j2012bl ). by 
considering the PC-current alone, also have predicted above behavior in the absence of strong 
rotation effect. Howe v er by considering effect of pulsar rotation on the field line constrained 
plasma alone, iBCWl (119911) have predicted that the lea d ing si d e inte n sity dominate ove r 



the trailing side and l ater 



Thomas fc Gangadharal (bp07h . bvksl (b()08h . byks et all feoioh 



Thomas et al. ( 2010 ). KG (j2012a ). Wang et al. ( 2012 ) have arrived at the same conclusion 



The "partial cone" pulsars, which show either missing or much suppressed side in their conal 
double component profiles (ILyne fc Manchester! Il988l ) can be result due to the rotation and 
PC-current perturbation (see Figures [131 and Ua 



Further the rotation and PC-current perturbation introduce an asymmetry into the 
width of components between the leading and trailing sides in such a way that the leading side 
components become broader than the corres ponding trailing ones, and it has an observational 
support too (lAhmadi fc Gangadharal 120021 ) . This is due to an induced asymmetry in the 
steepness of radius of curvature with respect to the r otation phas e w herein the tr ailing side 
p becomes more steeper than that on the leading side. IKGl (|2012aJ) and IKGl (l2012bl ) have also 
predicted the same behavior by considering the effects of rotation and PC-current separately. 

The rotation and PC-current significantly introduce an asymmetry into the strengths 
of opposite polarities of the circular polarization. When they are considered separately, they 
selectively enhance either the negative or positive polarities. Further the two effects cause 
the emission pattern to rotate in (6, 0)— plane in the opposite directions for positive a and in 
same direction for negative a. Hence it results in smaller rotation of the emission pattern in 
(6, cf))— plane for the positive a case than that for the negative a. However due to the rotation 
of the emission pattern, the symmetric type circular polarization becomes evident in addition 
to the more common antisymmetric type with an usual asymmetry between th e opposite 
polarities in presence of nonuniform distribution of sources (see 

alsoHEljSbl)- 

Due to 

perturbations we find the sign reversal of circular polarization from negative to positive near 
the center of pulse profile can be associated with either the increasing or decreasing PPA 
swing, and similarly the vice versa (see Fi gures [TTI and IT21 . and t hese deductions are found 
to be in accordance with lHan et all (Il998[ ) and lYou &: Hani (120061 ) observational findings. 
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Our simulations (see Figures [13] and [T4|) also confirm IHan et all (119981 ) and lYou fc Han 



(120061 ) findings that the negative polarity of the circular polarization is associated with 
increasing PPA and vice versa for positive polarity on both the sides of conal-double pulsars 
profiles. This is possible due to combined rotation and PC-current perturbation under 
specific conditions that when th e modulation is more effective in the polar direction than 
in azimuthal direction. However iKGl (j2012al ) have claimed that such correlation can exist 
provided in addition to the effective modulation in the polar direction, their location must 
be asymmetric in a given conal ring centered on the magnetic axis. 

Some pulsars particularly millisec ond pulsars show the polarization angle behavior that 
deviates from the standard 'S' curve (IXilouris et al.l 119981 ). We find that the 'kinky' type 
distortions in PPA profile are due to the superposition of modulated incoherent emissions 
over the regions which are larger than the wavel ength of rad i ation in p resence of s trong 
rotation and PC-current perturbations. Note that IGangadharal (120101) andlKGl (J2012a() also 
have obtained the similar results. However, the results of IGangadharal (120101 ) were in the 
absence of any perturbation and h ence the PPA distortions from the standard 'S' curve are 
smaller. iMitra fc Seiradakisl (120041 ) have speculated that if the radio emission has a varying 
emission height across the pulse profile then PPA swing will be nonuniform. Also, in the 
literature the kin ky type origin has been attributed to multi polar magnetic field in the radio 



emission region jMitra et a 



(IRamachandran fc Kramer 



2003 . 



2000) and to the return currents in the pulsar magnetosphere 



Note that we have performed modeling of coherent pulsar radio emission in terms of 
single particle curvature radiation. Although to the first order it is not a bad assumption, 
in reality the efficiency of coherence process may depend on many factors and hence there 
may be an induced asymmetries in the pulsar emission and polarization between the leading 
and trailing sides. For example the coherent emitters may be brighter on one side of the 
profile than on the other side depending upon the factors influencing the coherence process. 
And this can act in the same direction or in opposite direction to the processes that we have 
considered in this paper. 

Note that we have performed modeling of coherent pulsar radio emission in terms of 
single particle curvature radiation. Although to the first order it is not a bad assumption, 
in reality the efficiency of coherence process may depend on many factors and hence there 
may be an induced asymmetries in the pulsar emission and polarization between the leading 
and trailing sides. For example the coherent emitters may be brighter on one side of the 
profile than on the other side depending upon the factors influencing the coherence process. 
And this can act in the same direction or in opposite direction to the processes that we have 
considered in this paper. 
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5. CONCLUSION 

We believe our model is much more realistic for pulsar emission and polarization than 
those proposed before, as for the first time it simultaneously takes into account of the per- 
turbation caused jointly by the rotation and the PC-current. In addition it considers the 
detailed study on nonuniform source distribution (modulation) and the influence of viewing 
geometry on pulsar emission. Hence we believe one can explain most of the observed emis- 
sion and polarization properties of pulsars within the frame work of curvature radiation. On 
the basis of our simulation of pulse profiles we draw the following conclusions: 

1. The effect of rotation and PC-current perturbation in the presence of nonuniform 
source distribution (modulation) along with viewing geometry might be responsible 
for the most of the observed diverse behavior of polarization properties of pulsar radio 
emission. 

2. Because of the perturbations there arises an asymmetry in the phase shift of intensity 
components and PPA inflection point for ±cr. Further we notice that the phase delay of 
the PPA inflection point with respect to that of the central component (core) becomes 
larger for negative a than that for the positive a. 

3. The leading side components can become either stronger or weaker than the corre- 
sponding trailing side components of a given cone. This is due to the induced asym- 
metry in the curvature of source trajectory and the sight line encountered asymmetry 
in the modulation strength between the two sides. 

4. Both the "antisymmetric" and "symmetric" —types circular polarization are possible 
within the frame work of curvature radiation when the perturbation and modulation 
are operative. 

5. In the presence of perturbation the sign reversal of the "antisymmetric" —type circu- 
lar polarization as well as the sign of the "symmetric" —type circular polarization do 
not show the correlation with the sign of the PPA swing in the case of central core 
dominated pulsars. 

6. The 'kinky' type distortions in the PPA swing could be due to the incoherent addition 
of modulated emissions in presence of strong perturbations. 
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Fig. 1. — Geometry of curvature radiation in a rotating PC-current-perturbed magnetic 
dipole with the Cartesian coordinate system— XY Z as an inertial observer's frame whose 
origin is the neutron star center at O. The magnetic axis m is inclined to the rotation axis fl 
by an angle a, and the rotation phase 0' of the magnetic axis is measured from the fiducial 
XZ— plane. The thick solid line curves represent the PC-current-perturbed field lines, and 
the parameters a = 30°, 0' = 30°, field line constant r e = lr LC , pulsar rotation period 
P = 1 s, a current scale factor q — 1, magnetic azimuth 0« from 0° to 360° in steps of 45°, 
magnetic colatitude &i = 0°, and emission altitude r from to 0.8 r LC are used to sketch 
them. The unit vectors b and e represent the direction of perturbed field line tangent and 
pulsar rotation, respectively. The source net velocity -0 represent the direction of beamed 
radiation. 
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Fig. 2. — Emission point coordinates 9q and 0o as functions of rotation phase <$' in the cases of 
nonrotating dipole (thin solid line curves), nonrotating PC-current-perturbed dipole (dotted 
line curves), rotating dipole (dashed line curves), and rotating PC-current-perturbed dipole 
(thick solid line curves). The arrows in the 4>q panels (line styled similar to 0o) indicate the 
antisymmetric point of O in the respective cases. Here we used r„ = 0.1, P = 1 s, 7 = 400, 
and q — 1. 
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Fig. 3. — Radius of curvature p of source trajectory, and position angle ip of the radiation 
electric field as functions of <$' . The lines representation is the same as in Figure [2j The 
dot-dashed line curves represent the parameters obtained by simple addition of the two 
perturbations due to the rotation and the PC-current. The thin dot-dashed line curve 
in the ip panels represents the analytical perturbation theory fjHibschman &: Arond l200ll ) 
whereas the thick dot-dashed line curve represents our numerical result. The arrows in the 
i/j panels (styled same as if)) represent the inflection point of if). For simulation we used the 
parameters which are the same as in Figure EJ 
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Fig. 4. — The pattern of the circular polarization for the emissions from the beaming region 
with the uniform distribution of sources in the cases: the nonrotating dipole— panels (a) and 
(a'); rotating dipole— panels (6) and (6'); nonrotating PC-current-perturbed dipole— panels 
(c) and (c'); and rotating perturbed dipole— panels (d) and (cf). In each panel, emissions 
are normalized with the corresponding maximum of the total intensity. Here, we used the 
parameters = 0°, v = 600 MHz and the rest are the same as in Figure El 
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Fig. 5. — Simulations showing the emission with uniform distribution of sources before and 
after considering the perturbations due to the rotation and PC-current (the lines style is 
the same as in Figure [2]). In / panels, the emissions are normalized with the corresponding 
maximum of the total intensity after combining the the two perturbations, whereas in L and 
V panels, they are normalized with the corresponding maximum of I in the respective cases. 
The arrows in the i/j panels (styled the same as ip) indicate the PPA inflection points. Used 
v = 600 MHz and the rest are the same as in Figure El 
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Fig. 6. — Simulated pulse profiles with modulation in the azimuthal direction. The total 
intensity / is normalized with the maximum of / of the combined case of rotation and 
PC-current, whereas L and V are normalized with the corresponding maximum of / in the 
respective cases. We used r n = 0.05, fg — 1, <f)p — 0°, and = 0.1. The rest parameters 
are the same as in Figure [5j 




7. — Same as Figure [6] except r n = 0.1, a = 5° and two types of modulation: narrower 
= 0.1) and broader (o^ = 0.3) Gaussians are considered. 
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Fig. 8. — Same as Figure |6] but with two symmetrically located Gaussian modulations with 
respect to magnetic meridian plane and only the combined case of rotation and PC-current 
perturbations is presented. In the upper panels: p— dotted line curves, /—solid line curves; 
in the intensity panels: I— solid line curves, L— dashed line curves, V— dotted line curves. 
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Fig. 9. — Same as Figure [6] but with the nonuniform distribution of sources in the polar 
direction and = 1, Op = 2°, and oq = 0.01. 
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Fig. 10. — Same as Figure [9] with the combined effect of rotation and PC-current (the line 
style is the same as in Figure E]) except with the parameters 6 P = 3°. 6 and a e = 0.002. 
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Fig. 11. — Simulated pulse profiles with the nonuniform distribution of sources in both the 
polar and azimuthal directions, after considering the perturbations by the rotation and PC- 
current. The line style is the same as in Figure El We chose a = 5°, r n = 0.1, fo = 1, 
9 p = 3°. 4, (f)p = 0°, and the rest parameters are the same as in Figure 
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o- e =0.01, o- =O.l 0-^=0.001, 0-0=0.5 
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Fig. 12. — Same as Figure ITT1 except with a = 



-5° and <p P = 180°. 
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Fig. 13. — Same as Figure ITT1 but with two Gaussian modulations symmetrically located in a 
given conal ring at various azimuth. For simulation we used <jg = 0.003, = 0.1 and the rest 
parameters are the same as in Figure [TlJ The dashed line curves in the first two columns 
of the upper panels represent an amplified / to unity to identify with the corresponding 
intensity component. 
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Fig. 14. — Same as Figure [13] except with a = —5° and two Gaussian modulations considered 
in the first two column panels whereas four Gaussian modulations in the third and last 
column panels with different modulation parameters. 



